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1. SUMMARY 

A computer program, named INSL, is written for the calculation 
of the source location of impulsive noise. The main program is listed 
in Appendix 1. The algorithm is basically solving the triangulation 
equation: 

r^ + A.,(ifi)r 3 + A 2 (i^> + A 1 (ifi)r + A q ()|j) = 0 (1) 

where 

AqOJi) - [E - (T - A) U 0 F] 2 - E (T - A ) 2 C 2 

= 2 [2A - (T 1 - A) U q B] [E - (T - A) U q F] 

-2A(T - A ) 2 c 2 

A 2 W = [2A - (T - A) Uq B] 2 + 2[E - (T - A) U q F] 2 
- (T - A ) 2 c 2 

w 

A 3 (4») - 2 [2A - (T - A) UqB ] 

A = (-y^ sin ip + g sin a g cos i/; + cos ip) cos 
+ (g cos a g + Zj) sin g Q 

E = x.^ 2 + y^ 2 + z^ 2 + g (g + 2x 1 sin a g + 2z^ cos a g ) 

F = -g sin a - x. 

s 1 

B = -cos Bq cos i p 

C « 1052 + 1.143 TEM 

3 0 = 0 - a cos ip 

6 = a - a. cos $ - b. sin t p 
o Is Is 
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and 

r,\|) - Coordinates of source in ft and radian, respectively 

T = Time between acoustical signature and iJj = 0 index 

Uq = Tunnel velocity in ft/sec 

ti) = Rotor rotational velocity in Rad/ sec 

g = Shaft length above pivot in ft 

a = Shaft angle of attack in radian 
s 

a Q = Coning angle in radian 

a^ s = Longitudinal flapping coefficient 

b^ g - Lateral flapping coefficient 

TEM = Tunnel temperature in deg F 

x, y, z = Coordinates of microphone in ft 

The coordinates and symbols are illustrated in Figure 1. For detail 
formulation and discussion of the triangulation techniques, one should 
be referred to NASA CR 151996. 

For each azimuthal location a library subroutine is used to solve 
this fourth order equation, and all four roots are obtained. The algorithm 
picks out the meaningful roots, i.e., those which are Inside the rotor 
radius, then these roots are double-checked to ensure their accuracy. The 
procedures are then repeated for a different value of 1 p. Variable step 
sizes of ip are used to accelerate the computations. 

The computer program is written in FORTRAN for CDC 7600. The 
inputs are rotor operating conditions and the time intervals (T) between 
rotor l/rev index and impulsive noises as measured by different micro- 
phones. The outputs are the possible noise source locations in terms 
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of rotor radial and azimuthal coordinates. Typical computer time 
for a run of six microphone measurements is 1.5 sec, and the cost is 
about 12 cents for the CDC 7600 at NASA Ames Research Center. 
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ALGORITHM 


The block diagram is shown in Figure 2. The flow chart is shown 
in Figure 3. The coefficients of equation (1) are first calculated 
at an azimuthal angle ip based on the input data. A library subroutine 
ZPOLR is then called to solve equation (1) for radial location r. 

Four roots are obtained. They may be complex, real or mixed. Only 
positive real roots whose value are within rotor radius are meaning- 
ful, and those values are automatically selected. The value of <|i is 
then decreased by one step, and the calculations proceed. 

Two step sizes of are used. A coarse step size is used initially 
until a meaningful r is obtained. The azimuth, Tp, is then backed up one 
step size, and calculation proceeds with a fine step size. After a 
succession of calculation, there is a value of i(j beyond which no meaning- 
ful r is obtained. The calculations with fine step size persist beyond 
this point for several steps which is controlled as an input. If there 
is still no meaningful r obtained, the calculation proceeds with coarse 
steps until the range of ^ is reached, or until a meaningful r is obtained 
again. In the former case, the calculations re-initiate for the next 
microphone measurement. In the latter case, the step size is changed 
again in the manner previously described. 

The scheme of variable step size can save a large amount of computer 
time. Typically, the course step is 2 degrees and the fine step is 
1 degree for the first calculation. Depending on microphone locations 
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and other parameters, source lines extend either from the outer edge 
of the disk to the rotor hub or terminate at the middle. When calcu- 
lated results show source lines which terminate at the middle of the 
rotor disk, it may be real or may be due to the improper step size in 
ip. A smaller step size should then be used to recalculate for the 
source line within this rather small range of azimuth. For example, 
some calculations with 0.1 degree step size within an azimuthal range 
of 2° have been performed to obtain results shown in NASA CR 151996. If 
the source lines really terminate in the middle of the rotor disc, 
subsequent reduction of step size will not/change this result. 

The range of azimuth for calculation should be noted. For two 
bladed rotors, there are two acoustical impulses for each rotor period 
which is marked by l/rev index. The source location can be either 
before ip = 0 (1/rev index) or after it. The azimuth range of calcu- 
lation should cover both positive and negative azimuth. Either of the 
two acoustical impulses in a period can be chosen for calculation, 
although only one yields meaningful results. This situation must be 
resolved by physical arguments based on the characteristics of the 


noise source. 
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3. DATA CARDS 


A total of eleven (H) input data cards are used. 


Card Number Format Remarks 

1 2A1Q 20 alphameric characters used for 

labeling of output 

2 F10.3 Rotor angular velocity in rad/sec 

F10.4 Rotor shaft tilt angle in rad/sec, 

positive shaft tilt aft 

F10.4 Tip path plane angle in radians, 

positive leading edge up 

F10.1 Tunnel temperature in °F 

F10. 3 Tunnel velocity in ft/sec 

110 Number of microphones. Any number 

up to 6 can be specified 

3 6F10.4 X-coordinates of microphones, in ft 

4 . 6F10.4 Y-coordinates of microphones, in ft 

5 6F10.4 Z-coordinates of microphones, in ft 

6 6F10. 4 Time (in sec) between the 1/rev index 

and impulsive noise as measured by 
microphones. The values should be 
consistent with the sequence of micro- 
phones specified in Cards 3,4, and 5 

F10.4 
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F10.4 


Coarse step size of ij> in radius 
Fine step size of ij> in radios 


I 


7 

8 F1Q.4 Criterion of selecting real roots 

F10.4 Criterion of distinguishing double 

roots 

F10.4 Criterion in the double-check 

110 Number of calculation with fine increment 

after a succession of obtaining meaning- 
ful r 

9 F10.4 Largest value of azimuth ^ in radians 

F10.4 Smallest value of azimuth ijj in radians 

10 F10.2 Length of shaft in feet 

F1Q.4 Coning angle in radians 

11 F10.4 Radius of rotor in feet 
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4. EXAMPLE 

As an example, program XNSL is applied to a data point of high 
speed rotor noise measured in the 40- by 80-Foot Wind Tunnel at the 
NASA Ames Research Center. 

The rotor, test conditions, and input criteria are tabulated 
below: 


Rotor and Test Conditions 

Number of blades 

Blade precone angle 

Rotor radius 

Blade chord 

Blade twist 

Airfoil 

Hub articulation 

Shaft height above the pivot 
(See Fig. 1) 

Tip speed ratio, V/flR 

Rotor rotational speed 

Rotation tip Mach number 

Tip speed 

Shaft angle of attack 
Tip path plane 
Tunnel temperature 


2 

2 deg (0.0349 rad.) 

24 ft (7.315 m) 

1.75 ft (0.5334 m) 

-8 deg 

NACA 0012 

teetering 

6.94 ft (2.115 m) 

0.3 

327 RPM 
0.747 

821.8 ft/s (250.5 m/s) 

0 deg 
0 deg 

70 F (21 C) 

247.1 fps (75.32 m/sec, 146.4 knots) 


Tunnel speed 
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Microphones 

Number of microphones 

Microphone locations 


Time between 1/rev index 
and acoustical index 
(in sec) 

Input Criteria 

Coarse step 

Fine step 

Criterion of real root 
Criterion of double root 
Double-check criterion 
Largest ifi 
Smallest ip 
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X 


Y 


Z 

(ft) 

On) 

(ft) 

On) 

(ft) 

On) 

68.24 

20.80 

-20 

6.096 

8.06 

2.457 

6.99 

2.131 

5. 

1.524 

1.25 

.3810 

68.58 

20.90 

0 

0 

8.07 

2.460 

13.84 

-4.215 

10.92 

3.328 

-1.08 

.3292 

68.83 

20.98 

20. 

6.096 

8.08 

2.463 

12.41 

3.783 

12.83 

i 3.911 

7. 

2.134 

fit 

.1390, 

.0565 

, 1300, 

.0235, 

.1255, 

.0610 


2 deg 

(.0349 

radian) 



/ 

1 deg 

(.0175 

radian) 




1 

.5 

.001 

6.94 rad 
-3.1416 rad 


The input data deck to program INSL is shown in Appendix 2. 

The printout of the calculation is shown in Appendix 3. The first 
portion is the printout of input data. It then proceeds to print the 
results of the calculations. The first two columns are the coordinates 
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of source line in terms of their radial location r and azimuth location i p. 
DIS is the distance between source and microphone. TIME is the time for 
sound to travel from source to microphone, including effects of tunnel 
wind. CTIME is the calculated time between l/rev and acoustical impulse 
and should be consistent with input time (T) . A plot of source lines in 
this sample calculation is shown in Figure 4. 
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